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Abstract 

We introduce a new approach to analyse the global structure of 
electronic states in quasi- ID models in terms of the dynamics of a 
system of parametric oscillators with time-dependent stochastic cou- 
plings. We thus extend to quasi-lD models the method previously 
applied to ID disordered models. Using this approach, we show that 
a "delocalisation transition" can occur in quasi- ID models with weak 
disorder with long-range correlations. 

Pacs numbers: 73.20. Jc, 72.15.Rn, 05.40.-a 

1 Introduction 

In recent years there has been a steadily increasing interest in disordered 
models with long-range correlated disorder. The interest was initially spurred 
by the discovery that specific long-range correlations can produce a kind of 
"delocalisation transition" even in one-dimensional (ID) models in which the 
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electronic states are typically exponentially localised jl]. The first numer- 
ical results were later confirmed by analytical studies, which identified the 
relationship existing between localisation length and pair correlators of the 
random potential and showed how to create mobility edges in strictly ID dis- 
crete models [2 . These analytical predictions were experimentally verified 
by considering the transmission of microwaves in a single-mode waveguide 
with a random array of correlated scatterers jB]. Later on, the results ob- 
tained for discrete lattices [2J were extended to ID continuous models llj and 
applied to related problems such as the propagation of waves in waveguides 
with random surface scattering, and to specific quasi-lD models with bulk 
scattering |31IHIIS1- 

One of the main tools for the analysis of such models is based on the 
analogy between localisation phenomena in disordered systems and the dy- 
namics of stochastic oscillators [3 El- This approach was originally applied 
to strictly ID models; subsequent research, however, has begun to explore 
models of higher dimensionality, with quasi-lD models providing the first 
natural extension towards realistic disordered systems. A rigorous math- 
ematical treatment of discrete quasi-lD models was given in [9 , but the 
analysis was centred on disorder without spatial correlations. 

This paper serves two main purposes: in the first place, we extend to 
systems of higher dimensionality the dynamical approach which was used 
successfully for strictly ID continuous models (see, e.g., HTi] and references 
therein). We show that the Lyapunov exponents which govern the expo- 
nential divergence of initially nearby trajectories in a classical systems of 
stochastic oscillators are the same exponents which appear in the transfer- 
matrix approach for quasi-lD disordered models. This allows us to establish 
a rigorous correspondence between the quantum phenomenon of localisation 
in quasi- ID disordered models and the orbit instability in classical systems 
of parametric oscillators with noisy couplings. 

The second objective of this work is to use this analogy to study the effects 
of long-range correlations of the disorder on the localisation of the electronic 
states in quasi- ID models. Our main result is that, for weak disorder, specific 
long-range correlations can make all Lyapunov exponents vanish (within the 
second-order approximation), thereby suppressing orbit instability on the one 
hand and producing a "delocalisation transition" on the other. 

The paper is organised as follows: in SecElwe discuss the correspondence 
of quasi-lD disordered models with a set of coupled parametric oscillators. In 
Sec. ISlwe show how the evolution of this dynamical system can be analysed. 
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The general results thus obtained are then applied in Sec. El to the specific 
case in which the random potential depends only on the longitudinal coor- 
dinate. In Sec. El we apply the dynamical approach to the case of a generic 
weak disorder. We determine an expression for the sum of the positive Lya- 
punov exponents of the quasi-lD model and we use this result to discuss the 
delocalisation effects that are produced by specific long-range correlations of 
the random potential. The conclusions are then outlined in Sec. El 

2 Classical representation of the quasi- ID model 
2.1 The ID case 

Before considering quasi- ID models, we summarise shortly the main results 
for the strictly ID case. Since this work is focused on quasi-lD models, 
we shall be brief; the interested reader can find more details in [71 (HI ■ 
The correspondence between Anderson localisation in ID models with weak 
disorder and energetic instability of oscillators with a frequency perturbed 
by a noise is a straightforward consequence of the mathematical analogy 
between the Schrodinger equation 



with positive energy E and the dynamical equation of a stochastic oscillator 



(here and in the following we will use energy units such that h'^/2m = 1). 
In fact, Eq. (0) can be easily transformed into Eq. Q by interpreting the 
spatial coordinate x as the time t and the wavefunction amplitude i(j as the 
coordinate of an oscillator. In the first equation the function U plays the role 
of a random potential while in the second it represents a noise. The noise is 
white or coloured depending on whether the disorder is spatially correlated 
or not. 

The mathematical identity of Eqs. (PJ) and (0) allows one to study the 
global structure of the quantum eigenstates of the disordered model (|H) by 
analysing the dynamics of the corresponding classical oscillator The 
dynamical equation (0) gives the time evolution of the oscillator coordinate 
q{t) once the initial position g(0) and velocity g(0) have been specified; use of 




(1) 




(2) 
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analogous boundary conditions to solve the Schrodinger equation ((T)) leads 
to the solution iIj{x) which is obtained with the standard transfer matrix 
approach. Studying the disordered model (Q) in terms of the dynamics of 
the random oscillator therefore, is equivalent to using transfer matrix 
methods. 

Comparing the solutions of Eqs. (Q) and Q one finds that spatially ex- 
tended states correspond to bounded oscillator orbits, while localised states 
have their counterpart in unbounded trajectories. As a consequence, the phe- 
nomenon of Anderson localisation corresponds to energetic instability of the 
parametric oscillator, with the inverse localisation length being equal to the 
Lyapunov exponent of the stochastic oscillator, i.e., the rate of exponential 
divergence for initially nearby orbits. 

The case of weak noise/disorder can be studied using perturbative tech- 
niques which were originally devised for the study of stochastic systems. 
These methods allow one to obtain the rate of energy growth for the os- 
cillator and therefore the inverse localisation length for the disordered 
model ([[]) [HIE]- The second-order expression for the inverse localisation 
length is 

A = — / {U{t)U{t + T)) COS {2uT)dT (3) 


which shows that the Lyapunov exponent is proportional to the power spec- 
trum of the disorder, i.e., to the cosine Fourier transform of the two-point 
correlator of the random potential (here and in what follows we use the sym- 
bol (...) to denote the average over different disorder realisations). This 
entails that the energetic instability of the oscillator Q can be suppressed 
if the unperturbed oscillator frequency, multiplied by a factor two, lies in a 
frequency interval where the power spectrum of the disorder vanishes. Corre- 
spondingly, delocalised states arise for the energy values for which the inverse 
localisation length (jH)) vanishes. 

2.2 The quasi-lD model 

We analyse the phenomenon of electronic localisation in a semi-infinite strip. 
The present method can be applied to bars as well; we focus on strips to 
avoid unnecessary complications in the mathematical formulae. We consider 
the strip D in the x — y plane 

D = {{x,y):0<x; 0<y<L}. (4) 
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The Schrodinger equation for a quantum particle ( "electron" ) in the strip is 




^ + ^ ^) + y)^(^^ y) = Ei){x, y). (5) 



The function U (x, y) is the random potential felt by the electrons; the po- 
tential can exhibit spatial correlations. The dimensionless parameter e is 
introduced to keep track of the order of perturbative expansions and can be 
set equal to 1 in the final formulae. We will focus on the case of weak disor- 
der, i.e., on the case e <^1. For the model to be completely defined, one must 
provide the statistical properties of the random potential; for weak disorder, 
it is enough to specify the first two moments of U{x,y). We will assume 
that the potential has zero average, {U{x,y)) = 0, and that the two-point 
correlator is a known function 



provided that it is finite; when the variance diverges, like in the case of white 
noise, must be interpreted as a parameter which measures the strength of 
the disorder. We will consider random potentials which are translationally 
invariant in the mean; this is why we assume that the binary correlator 
depends only on the difference of the coordinates of the two points {xi,yi) 
and {X2,y2)- 

2.3 Representation of the Schrodinger equation in the 
partially reciprocal space 

If one wishes to extend to quasi-lD models the dynamical approach described 
in subsection 12 . 1 I for the ID case, it is quite natural to identify the longitudinal 
coordinate x with the time t of the corresponding dynamical system. One 
can then Fourier-transform the Schrodinger equation (0) in the transversal 
direction y. In the following we consider the real solutions of the stationary 
Schrodinger equation (j3J; this choice is not restrictive because one can find 
a complete set of real eigenfunctions of the Hamiltonian. It is possible to 
expand any real function of y defined over the interval [0 : L] in a cosine 



{U{xi,yi)U{x2,y2)) = o-\(a;i 



(6) 



where cr^ 



represents the variance of the potential 



= {U{x,y)U{x,y)) 
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basis. In this way one avoids deahng with complex Fourier components, 
which makes easier the successive identification of the Fourier components 
with oscillator coordinates discussed below. Expanding the wavefunction and 
the potential one obtains 



oo 

Tiny 



L 

. (7) 
7my^ 



and 



i'{x,y) = ^o{x) + 2^ 'ipn{x) cos (^^^^ = ^ '(pn{x) cos 

n=l n=— oo 

oo oo 

U{x,y) = Uoix) + 2 J2Un{x) cos (^^^ = ^[/„(x)cos( 

n=l n=— oo 

where the Fourier components of the expansion are defined by 

1 fTTny\ 
^n{x) = -j dy ^{x,y) COS [-J- j (8) 

Unix) = - J dy U{x,y) COS {-j-y (9) 

In the following we will often refer to the Fourier components of the wave- 
functions as Fourier modes or channels. Note that the Fourier components 
with opposite values of the index n are identical: ipn^x) = and 
Un{x) = U^n{x). The functions defined by the Fourier cosine expansion ((Tj) 
coincide with the original functions in the domain outside of this domain, 
they are symmetric and 2L-periodic functions of the transversal coordinate 
y, i.e., 

^Ij{x,y) = ^lj{x,-y) and ip{x,y) = ipix.y + 2L). 

The same relations also hold for U{x, y). Expanding in a cosine series, there- 
fore, is equivalent to consider ?/-symmetric functions ip{x,y) and U{x,y) on 
the doubled strip [0 : oo] x [—L : L] and impose periodic boundary conditions 
along the transversal direction (i.e., to roll the doubled strip into a cylinder). 
Boundary conditions can be chosen with considerable latitude because they 
do not infiuence the structure of the electronic state if the strip is sufficiently 
wide. By taking the Fourier cosine transform, one can write Eq. (0) in the 
form 

^ ^"^a;) + E- (^) V^„(x) = J2 eUn-k{x)Mx). (10) 



dx"^ 
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2.4 Elliptic and hyperbolic Fourier components 

Before giving a dynamical interpretation for Eq. (fTUI) . it is useful to intro- 
duce a distinction between "elliptic" and "hyperbolic" Fourier components 
of the wavefunction. Adopting the terminology of , we define the Fourier 
component ipniu) with wave number n to be 



eUiptic if E — 
hyperbolic if E — 




(we ignore the marginal case of the "parabolic" components with E = 
(nn/Lf). 

The crucial difference between elliptic and hyperbolic components is that 
the latter decay exponentially in the longitudinal direction: from a physical 
point of view, they are evanescent modes (or closed channels). Hence they can 
be neglected for large values of the longitudinal coordinate x. The irrelevance 
of the hyperbolic modes can be justified also with a different argument. As 
can be seen from the Fourier expansions (|7j), the Fourier components with 
large values of n describe the behaviour of the wavefunction on small spatial 
scales ^ 

Sy 

n 

in the transversal direction. If we suppose that our continuous model has an 
underlying lattice structure, however, we need not consider the behaviour of 
the wavefunction over spatial scales Sy < a, where a is the lattice constant. 
This implies that, if the short-scale spatial structure of the wavefunction is 
neglected, the only Fourier components that must be taken into account are 
the ones with wavenumber 

_ L 

n < n —. (11) 

a 

Now, let us restrict our analysis of localisation to the case of high-energy 
electrons, where high-energy means 

In this case, one can see that, due to the cutoff condition (fTTj). the only Fourier 
components which are relevant for the analysis are the elliptic ones. Hence, 
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we need not bother with hyperbohc components if we restrict our attention 
to the regime of sufficiently high energies. This back-of-the-envelope criterion 
suggests that our approximation may fail when the electron energy does not 
satisfy condition (fT^ . 

After eliminating the hyperbolic components, the Schrodinger equation ITUI) 
takes the form 



x) + 



N 



tpn{x) = ^ eUn-kix)lpk{x) 



(13) 



-N 



where the indices n and k are restricted to the elliptic modes (or conducting 
channels), i.e., the indices take the values — A^, — + 1, . . . , A^ — 1, A^ with A^ 
being the integer part [. . .] of the ratio L^/E/tt, 



N 



TT 



2.5 Dynamical equations for the quasi-lD model 

In order to give a dynamical interpretation of the Schrodinger equation (jl3|) 
we replace the longitudinal variable x with a time variable t. In addition, we 
define the frequencies 

and we introduce the notation 

qn{t) = 1pn{t) 

for the Fourier components of the wavefunction. This allows us to write 
Eq. (jl3|l in the form 

TV 

qn{t) + uolqnit) =e J2 Un-k{t)qk{t) (15) 

k=-N 

which is naturally interpreted as the dynamical equation of a system of clas- 
sical parametric oscillators with time- dependent stochastic couplings. This 
correspondence allows one to analyse the structure of the electronic states of 
the disordered model (0) in terms of the dynamics of the system (jl5j) . In fact. 
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the spatial behaviour of the ipn{x) Fourier component along the longitudinal 
direction is determined by the time evolution of the coordinate qn{t) of the 
corresponding oscillator. 

The system of 2N + 1 second-order differential equations (fT3j) can be 
transformed into a system of 4N + 2 first-order Hamiltonian equations by 
introducing the momenta p„ = g„ and the Hamiltonian 

N y 2 2 \ N N 

H{P,g)= E Y + Y^' ^ ^ Un-,{t)qnqk. (16) 

n=-N ^ ^ n=-Nk=-N 

Then one can cast the dynamical system (|TH)l in the form 

OH ^ 

Pn = -TT- = -^Iqn + ^ Un-k{t)qk 

k=-N (17) 

dH 

qn = = Pn- 

OPn 

Note that these Hamiltonian equations describe a system of 2N + 1 oscil- 
lators. However, not all oscillators are independent, because the symmetry 
of the Fourier components ipni^) = ip-n{x) implies that the Hamiltonian 
system is subject to the constraints g„ = g_„ which reduce to + 1 the 
number of degrees of freedom of the system. Making use of the constraints, 
the dynamical equations ()17|) can be written in the form 

N 

Pn = -^Iqn + £Wnk{t)qk 



fe=0 



qn = =Pv 



where n = 0,1, ... ,N and only independent oscillators are involved. Note 
that in Eq. (fTHj) we have introduced the short-hand notation 

Wn,k{t) = [Un+k{t) + f/„_fe(t)] (^1 - ^5k0^ . (19) 

A comment is in order here. A mapping of a quasi- ID disordered model 
of the form ^ unto a Hamiltonian system was made long ago by Hansel 
and Luciani fT^ llHj. In spite of several similarities, their approach differs 
from the present one in many aspects. In the first place, in fT^ IT^ the 
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quasi- ID model is morphed into a Hamiltonian system by discretising rather 
than Fourier-transforming the system in the transversal directions. In our 
approach, no discretisation is required. This implies that the channels Hansel 
and Luciani deal with are ID chains coupled linearly to each other, whereas in 
our scheme the channels are the Fourier modes of the wavefunction. Because 
of this conceptual difference, our channels are coupled only if the random 
potential is present, contrary to what happens in the Hansel and Luciani 
work. Similarly, no distinction between elliptic and hyperbolic channels arise 
in that case. In the second place, the Hamiltonian model of Hansel and 
Luciani has no constraints, while the system (fT^ is subject to the constraints 
g„ = g_„. This entails that the coupling term (jl9|) is not fully symmetric in 
the two indices n and k, unlike the coupling matrix considered in [T^ I13j. 

2.6 Wavefunction localisation versus oscillator insta- 
bility 

The mapping of Eq. (fT^ unto Eq. (fTHj) makes possible to study the behaviour 
of the wavefunctions of the disordered model (0) in terms of the time evo- 
lution of the dynamical system (jl8|) . From a physical point of view, the 
mathematical identity of the two problems translates into a correspondence 
between the phenomenon of Anderson localisation in quasi-lD disordered 
bars and the energetic instability of the system of random oscillators. The 
connection between the two phenomena can be understood and made quanti- 
tative by observing that the crucial properties of both classes of systems are 
defined in terms of their Lyapunov exponents. In the case of quasi-lD sys- 
tems, Lyapunov exponents arise within the framework of the transfer matrix 
approach (see, e.g., jT3] and references therein). The method divides the strip 
(or the bar) in layers; a transfer matrix is an operator that relates the values 
of the wavefunction and its derivative on one layer with the corresponding 
values on the contiguous layer. The Schrodinger equation is considered as 
an initial value problem and its solution is obtained in terms of a product 
of transfer matrices (in discrete models) or of spatial-ordered exponentials 
(for continuous models). Making use of Oseledec's theorem ^3], one can 
then define the Lyapunov exponents for the quasi- ID model in terms of the 
eigenvalues of the asymptotic product of transfer matrices. The localisation 
length is equal to the smallest positive Lyapunov exponent. 

For dynamical systems, the Lyapunov characteristic exponents can be in- 
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troduced by considering the exponential divergence of initially nearby orbits. 
For maps, this leads to the computation of the eigenvalues of limit product of 
matrices (see, e.g., [^) which give the evolution of vectors in tangent space 
and are the counterpart of the transfer matrices for the quasi-lD models. 
To sum up, Lyapunov exponents can be defined both for quasi-lD models 
and for dynamical systems; in particular, due to the mathematical identity 
of Eqs. (fT^ and (fTHj) . the (non- negative parts of the) two Lyapunov spectra 
are the same for both systems. 

From this point of view, quasi- ID models are similar to ID systems of 
the type described by Eq. ((H), whose single Lyapunov exponent (i.e., the 
inverse of the localisation length) is identical to the Lyapunov exponent of 
the dynamical counterpart (j2I) (i.e., the rate of exponential orbit divergence). 
However, in the case of quasi-lD systems, there is an important difference 
between the disordered model (0) and its dynamical analogue fITHjl : for the 
former, the most important Lyapunov exponent is the smallest one, which 
defines the inverse localisation length, while for the latter the dynamics is 
dominated by the largest Lyapunov exponent, which gives the rate of ex- 
ponential divergence along almost every direction of the phase space of the 
system. Therefore, in order to extract information on the spatial behaviour 
of the electronic states of model a thorough analysis of the time evolution 
of the dynamical system ()18p is needed. In particular, one cannot restrict 
the dynamical analysis to the determination of the largest Lyapunov expo- 
nent. Having thus highlighted the links between the quantum phenomenon 
of localisation and the dynamics of a system of parametric oscillators, we 
now turn our attention to the study of the latter problem. 

3 Dynamics of the parametric oscillators 

To analyse the dynamics of a system of oscillators, it is convenient to perform 
a canonical transformation and switch from the Cartesian coordinates (g„, p„) 
to the action-angle variables ( J„, 6'„) defined by the relations 




Pn = 



^/2uJnJn cos 9, 
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In terms of the new variables the Hamiltonian takes the form 



N N N I J J 

H= ^ uJnJn- ^ ^ £:f/„_fc(t)W^^sin6'„sin6'fc 

n=-N n=-Nk=-N ' ^n^^k 

and the dynamical equations (fTRjl become 

TV 



Jn = y^2eWn,k{t)\-^sin6kCO?,9n 

^ I — Jk — 

On = uJn-y^eWn,k{t)\ - — - — siu 6^^ siu 6'„ 

f— ^ V JnUJnUJk 

k=0 



(20) 



Following the general method described by Van Kampen one can 
replace the system of Langevin equations piU) with a deterministic Fokker- 
Planck equation for the probability distribution of the stochastic variables 
{Jn-iOn)- In fact, the dynamical equations (plj) are stochastic equations of 
the form 

M = F(°)(M)+£F(i)(M,t) (21) 

where F^'^^u) is a sure function of the vector u and eF^^^u, t) is a stochastic 
correction to the deterministic term. Following Van Kampen J7], one can 
associate to the stochastic differential equation (PT|) an ordinary differential 
equation for the function P{u, t) which represents the probability distribution 
for the random variable u. The associated differential equation has the form 

ij J 

(22) 

where m* is the solution at time t of the unperturbed differential equation 
ii = F^^\u) with the initial condition u{0) = u, and d{u~^) /d{u) is the 
Jacobian of the transformation u u~'^ . Note that, within the second-order 
approximation, i.e., neglecting terms of order 0(5^), only first- and second- 
order derivative appear in the right hand side (r.h.s.) of Eq. (f^^ . which 
therefore has the form of a Fokker-Planck equation. 
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In the present case, the vectors u, F^^\u), and F^^\u,t) have the n-th 
bi-component equal to 



where we have introduced the symbols 

N 



F^nliJ^d) = y]2W^„,fc(t)W^^sine,cos^^„ 

N I J 

Fi^l{J,0) = -y^Wn,k{t)\—^Mu^inen 

T~, V Jn^n'^k 

k=0 

The unperturbed flow m — in this specific case has the simple form 

^ P \ ( 1 



and as a consequence one has 

d (u-^) _ d{u) _^ d _ d 

d {u) d (W^) du~'^ dui 

Using these results, the general Fokker-Planck equation (|22j) takes the specific 
form 

BP ^ BP 

n=0 

N N r. „r^. 



n=0 fc=0 

N N 



n=0 fc=0 

N N 



n=0 A:=0 
N N 



UJT, t - 


-T))P{J,9,t) 


UT, t — 


T))P{J,9,t) 


UJT, t — 


T))P{J,9,t) 


UJT, t — 


r))PiJ,9,t) 



n=0 k=0 



(23) 
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Rather than trying to obtain the general solution of Eq. (|23|) -an exceedingly 
difficult task-, we aim at deriving the reduced probability distribution for the 
action variables 

PiJ,t) = I \[ddnP{J,9,t). 

Since the frequencies Un are different from zero, a glance at the dynamical 
equations (f^UI) shows that the angular variables are "fast" in comparison with 
the "slow" action variables. We can therefore assume that, after a sufficiently 
long time, the angle variables become uncorrelated random variables with a 
uniform distribution in the [0 : 27r] interval, i.e., that for large times the 
distribution P{J,9,t) reduces to the factorised form 

(24) 

Substituting the distribution in the Fokker-Planck equation (^Hj) and 
integrating over the angular variables, one obtains, after some algebra, the 
reduced Fokker-Planck equation 



BP " " 1 a 

n=0 k=0 " 



dP 

Dr.,{J) — {J,t) 
OJk 



(25) 



where the diffusion matrix Dnk{J) has diagonal elements defined as 

Dnn{J) = (^^) ^ dT {Wnn{t)Wnn{t-T)) COS 2UJ,,T 



+ 52 — dr {Wnk{t)Wnk{t-T)) [cos {uJn + UJk) T + COS {Un- UJk) r] 

(26) 



and off-diagonal elements (n 7^ k) equal to 

t 'Jn'Jk 



Dnk{J) = / dr {Wnk{t)Wkn{t-T)) [cOS (t^„ + Uj,) T - COS {uJn - UJk) t] . 

UJnUJk Jo 

(27) 

The Fokker-Planck equation is a rather general result, because it has 
been obtained for a generic random potential, with the only assumption that 
the disorder should be weak and translationally invariant in the mean. At 
ffist sight, one could find surprising that a localisation phenomenon can be 
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described with a diffusion equation. The apparent paradox can be explained 
by considering that the locahsed Fourier components of the electronic wave- 
function have oscillators with unbounded orbits as dynamical counterpart. 
This equivalence provides an intuitive reason why electronic localisation can 
be analysed with a diffusion equation for the action variables. 

We note that only second-order derivatives appear in the diffusion equa- 
tion (j25|) . and that the coefficients of its r.h.s. are completely determined 
by the two-point correlator of the random potential. Both facts are conse- 
quences of the decision to study the problem within the framework of the 
second-order approximation. Pushing the perturbative approach to higher 
orders, in fact, would make higher-order derivatives appear in the diffusion 
equation, together with coefficients which would depend on higher-order mo- 
ments of the potential. It may be appropriate to stress that, if in the present 
work the statistical properties of disorder are defined only up to the bi- 
nary correlator, this is not due to any Gaussian assumption but is only a 
consequence of the adopted second-order perturbative scheme, which makes 
unnecessary to specify the statistical features of the random potential beyond 
the second moments. 

Using the Fokker-Planck equation ()25|) for the reduced probability density 
P{J,t) allows one to obtain the dynamical equation for the averaged action 
variables 



In fact, differentiating with respect to time both members of the previous 
equation and using Eq. (^3)) to express the time derivative of P{J,t), one 
arrives at the matrix equation 




d 
dt 



I 



M 



(28) 



\Jn j 



\ 



Jn j 



where M is a (A^ + 1) x (A^ + 1) matrix with diagonal elements 



M 



-nn 
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and off-diagonal elements (n 7^ k) 

^2 POO 

Mnk = / dr {Wnkit)Wnkit-T)) [cOS (uj^ + t^fc) T + COS (u;„ - UJk) t] . 

ZOJn^k Jo 

An important consequence of Eq. (PH|) is that, except in very special cases 
(discussed in Sec. E}, the exponential rate of energy growth is the same for 
all oscillators. In fact, every action variable increases exponentially in time 
with a rate given by the largest eigenvalue of the M matrix. 



4 Longitudinal disorder 

As a first application of the previous results, we can consider the special case 
in which the random potential depends only on the longitudinal coordinate, 
i.e., U{x,y) = U{x). We will refer to this "longitudinal disorder" . 

This form of random potential has been considered in a different context, 
i.e., that of many-mode waveguides with a rough surface, where it has been 
christened as "stratified disorder" JH]- Here we recover the results of that 
paper using our more general formalism. 

In the special case in which U depends only on the longitudinal coordinate 
X, the matrix elements (fT^ take the simple form 

Wnk = U{t)5nk. (29) 

Substituting this expression in Eqs. (jSE)) and one obtains that the dif- 
fusion matrix becomes 

Dnk{.J) = 4A„J^5„fc 

where we have introduced the symbols 

A„ = — / dr {U{t)U{t + r)) cos 2c<;„r. (30) 



As a consequence, the Fokker-Planck equation reduces to 



N 



dt ' ^ dJn 



The solution of this Fokker-Planck equation corresponding to the initial dis- 
tribution 



N 



P(J,t = 0) = J](5(J„-J40)) 



?i=0 
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is a product of log-normal distributions 



N 



1 



[log(JjJ„(0))-4A„t]^ 

8Xnt 



p(j,t)=n 



exp 



(31) 



71=0 



Using expression ()29|) one can also compute the matrix elements of the 
M matrix which determines, via Eq. ()28|) . the time evolution of the averaged 
action variables. In the case of longitudinal disorder the M matrix takes the 
simple diagonal form 



Both this result and the factorised distribution (jHT|) imply that, in the case 
of longitudinal disorder, the parametric oscillators are decoupled. The re- 
sult could have been obtained a priori by noting that if the random poten- 
tial depends only on the longitudinal coordinate, Fourier transforming the 
Schrodinger equation (0) gives a system of independent equations for the 
Fourier components. Thus, the quasi- ID model (0) is effectively decomposed 
into + 1 strictly ID systems, in agreement with the results of [TH] . 

Besides being independent, the oscillators are also energetically unstable; 
in fact, solving Eq. one obtains 



The coefficients 4An are therefore the mean rates of exponential growth of the 
energies of the oscillators (since the energy of the n-th oscillator is propor- 
tional to the corresponding action variable. En = uOnJn)- Taking into account 
that the rate of exponential increase of the energy is four times the rate of 
exponential orbit divergence [8J, we are led to the conclusion that the coef- 
ficients flSUj) represent the Lyapunov spectrum of the dynamical system (fTH|l 
and, consequently, of the quasi-lD model (0). The localisation length, there- 
fore, is the equal to the inverse of the minimum Lyapunov exponent. Which 
of the Lyapunov exponents (j3(J|) is the smallest depends on the specific form 
of the Fourier transform of the binary correlator {U{t)U{t + r)) = (T^x(r). 

In the case of ^-correlated disorder, the Fourier transform of x(t) = 
(5(r), is simply the unity, xi^) = and the smallest Lyapunov exponent 
corresponds to the oscillator with largest unperturbed frequency, i.e., the 
oscillator with n = 0. Hence the inverse of the localisation length is 



1 




8E 
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When the binary correlation function is not a delta, however, the smallest 
Lyapunov exponent is not determined only by the largest frequency, but also 
by the behaviour of the power spectrum of the random potential 

/oo /»oo 
X(r) cos UTdr = — {U{t + r)f/(t)) cos tuTdr. (32) 
-oo ^ J —oo 

An important consequence is that the system can go through a delocalisa- 
tion transition if long-range correlations of the disorder make the Fourier 
transform of two-point correlator ()32|) vanish in a specific frequency interval. 
This phenomenon has already been analysed for ID systems [21 El ^] and 
for quasi-lD waveguides with stratified disorder [TH]. 

The main point is that, given a Lyapunov exponent Xn{uj) with any spe- 
cific dependence on the frequency u, it is possible to find a random po- 
tential U{t) that generates the pre-defined function A„(u;). Specifically, if 
the Lyapunov exponent A„(u;) is known, then the power spectrum x{uj) = 
Suj'^ Xn{uj) / a"^ is also defined. One can then determine a function j3{t) whose 
Fourier transform is a/^(cJ) via the inversion formula 

.duj 



Pit) = j ^^)e'^-' — . 

A random potential which produces the desired behaviour of \niyj) is then 
obtained by taking the convolution of the function (3{t) with a stochastic 
process ri{t) with zero mean and delta-shaped correlation function. In other 
words, one can consider the potential 



U{t) = / ds P{s)ri{t + s) 

J — oo 

where ri{t) is a white noise with 

{r]{t)) = and {r]{t)r]{t + r)) = S{t). 

Following this recipe, one can obtain for instance a random potential with 
the long-ranged correlation function 

x(t) = — (sin UiT — sin h'2T) (33) 
r 



which corresponds to the "window" power spectrum 



1 a Ui < CO < 1^2 

otherwise 



18 



In such a case, one has that the for every frequency ui^ that falls outside 
of the interval [z^i/2 : 1^2/2] the corresponding Lyapunov exponent van- 
ishes (at least within the limits of the second-order approximation considered 
here). By shifting the frequencies z/i and 1^2, one can therefore obtain a de- 
localisation transition as soon as the smallest Lyapunov exponent vanishes; 
this corresponds to the electronic wavefunction having one extended Fourier 
component. If more Lyapunov exponents vanish, the number of extended 
Fourier components increases and the delocalisation effect becomes more ro- 
bust; in the extreme case when all Lyapunov exponents vanish, the electronic 
wavefunction is not affected by the random potential. If the phenomenon is 
considered from the point of view of the dynamical system (|18|) one has that 
for every vanishing Lyapunov exponent there is an oscillator which becomes 
stable energy-wise. When all the Laypunov exponent are zero, all the oscilla- 
tors are stabilised, i.e., the exponential divergence of the orbits is suppressed 
and the dynamics ceases to be chaotic. The possibility of making selected 
Lyapunov exponents vanish entails that in solid-state models the transmis- 
sion properties can exhibit anomalous and unexpected features of selective 
transparency, as happens for waveguides. We refer the reader to ^H] for a 
detailed discussion of this phenomenon. 



5 Delocalisation transition in quasi- ID mod- 
els 

In the previous section we have considered the special case in which the ran- 
dom potential depends only on the longitudinal coordinate. The methods of 
Sec. 0] can be applied also to the a potential of the general form U = U{x, y), 
but in this case it is usually impossible to find the analytic solution P{J,t) 
of the Fokker-Planck equation as well as to solve with non-numerical 
methods the differential equation (j^Hj) for the average action variables. In 
fact, the dependency of the random potential on the transversal variable pro- 
duces coupling of the different oscillators, so that the coupling matrix (fT^ 
is no longer diagonal. As a consequence, neither the diffusion matrix in 
the Fokker-Planck equation (j^^ nor the evolution matrix Mnk in Eq- (j^Hj) 
are diagonal and this makes exceedingly difficult to solve analytically both 
equations. 

One can, however, use a different technique to determine the sum of the 
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positive Lyapunov exponents, or Kolmogorov entropy. Setting the latter 
equal to zero then gives a sufficient condition for delocalisation of the elec- 
tronic states of the quasi-lD model (0) and a sufficient and necessary con- 
dition for the suppression of orbit instability in the dynamical system ()18|) . 
In the following subsection we define the Kolmogorov entropy and we show 
how it can be computed in the second-order approximation. The reader un- 
interested in technical details may skip this subsection and the next, where 
a few special cases are analysed, and go to Subsec. 15.31 where delocalisation 
effects are discussed. 

5.1 The Kolmogorov entropy 

It is well known that in a Hamiltonian system with v degrees of freedom 
there are 2v Lyapunov exponents which, due to the symplectic structure of 
the dynamical equations, obey the symmetry relation 



with % = (see, for instance, Because of this relation, one 

has u non-negative Lyapunov exponents. For a deterministic Hamiltonian 
system, at least one of these exponents vanishes; in the present case, however, 
the presence of a noisy term in the Hamiltonian ensures that, under 
normal circumstances, the non-negative exponents are actually positive. The 
Lyapunov exponents are defined as the exponential rate of local divergence of 
initially nearby trajectories; this definition, however, can be used operatively 
only to compute the largest Lyapunov exponent. Lesser Lyapunov exponent 
can be determined using a technique devised by Benettin et al. [20| . The main 
idea is that the sum of the largest k exponents is equal to the exponential rate 
of increase in time of the volume of a parallelepiped spanned by k independent 
vectors ( with k <v). The result does not depend on the choice of the initial 
vectors. 

In the present case, we are interested in the sum of all the positive Lya- 
punov exponents; for this reason we consider the volume of the parallelepiped 
spanned hy v = N + 1 linearly independent vectors . . . , i^^\t). Such 

a volume can be expressed as the square root of a Gram determinant 



Aj — —\2u-i+l 



V^''\t) = ^\detG{t)\ 
where G is the (A^ + 1) x [N + 1) matrix with elements 



(34) 



(35) 
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The sum of the + 1 positive Lyapunov exponents can then be written as 

(36) 



i=0 



To apply this prescription to the dynamical system (fTHj) it is convenient 
to introduce the rescaled variables 

_ Vn 



with n = 0, . . . ,N, so that the unperturbed motion of the n-th oscillator 
reduces to a rotation in the {xn,XN+i+n) plane. 

In terms of the new variables, the dynamical equations pSj) take the form 



N 



Xr 



—U!nXN+l- 



(37) 



The time evolution of the system can be expressed in terms of the evolution 
operator U(t), defined by the relation 



(38) 



Let us introduce the (2A^ + 2) x (2A^ + 2) matrices A and B, defined in block 
form as 







( ^o(O) \ 




= U(t) 




\ a;2iv+i(t) j 




\ a;2jv+i(0) J 



A 



-n 
n 



and B(t) 



ri-i/2w(t)rj-i/2 





where W(t) is the (A^ + 1) x (A^ + 1) matrix whose elements are defined by 
Eq. (Uni) and ft is the (iV + 1) x (A^ + 1) diagonal matrix 



nk 



^n^nk- 



Writing the dynamical equations ()37|1 in matrix form and taking into account 
Eq. (jSHl), it is easy to see that the evolution operator is the solution of the 
matrix differential equation 



U = [A + £B(t)]U 
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with the initial condition U(0) = 1. Going to the interaction representation, 
one can write the evolution operator in the form 

U(t) = e^'U^t) 

where the first factor is the unperturbed evolution operator, 

( cos fit — sin Q,t 
\y sin Ot cos fit 

while the second factor obeys the equation 

U/ = eB/(t)U/ (39) 

with 

B,(t) = e-^*B(t)e^*. 
For our purposes it is useful to write the evolution operator in block form: 

This decomposition allows one to obtain from Eq. (jH^ the dynamical equa- 
tions for the left blocks of the evolution operator in interaction representation 



(40) 



where we have introduced the new symbols 



Wcc(t) = cos fit n-^/^ w(t) n-^/^ cos fit 

Wc,(t) = cosfit fi-i/2 W(t) fi-^/2 gj^Q^ 

W,c(t) = sin fit fi-^/2 W(t) fi-^/2 COS Ot 

W,,(t) = sin fit fi-^/2 W(t) fi-^/2 gjj^ 



The identity U7(0) = 1 implies that Eqs. ()40|) are to be solved with initial 
conditions 

U;"^(0) = 1 and Uf(0) = 0. (41) 
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We can now select as initial parallelepiped the (A^+ 1) -dimensional cube 
of edge A in the space of rescaled momenta. In other words, we consider the 
set of initial vectors 

(0) = A5,fc 

where the index i (which identifies the independent vectors) runs from to 
A^, while the index k (which labels the components of each vector) runs from 
to 2N + 1. With this choice of the initial vectors, and remembering that 
the evolved vectors can be written as ^'■'\t) = V{t)^^'\0), one obtains that 
the matrix (jH^j) can be expressed in terms of the left blocks of the evolution 
operator 



G 



A^ 



(where the symbol denotes the transpose of the matrix M). Note that 
the matrix G has the same form in the interaction representation and in the 
original representation. Inserting this matrix in formula fl34j) . one obtains 
the volume of the expanding parallelepiped; after substituting this result in 
expression (j36p one arrives at the conclusion that the sum of the positive 
Lyapunov exponents can be written as 



N 



lim — ( log det 

T-*oo 2T 



\jf^^{T)ijf\T) + ^(r)uf (T) 



lim — - 

T^oo 2T 



dt ( — log det 
\dt 



vf^' (t)Uf (t) + U^"^^ {t)Vf>{t) 



(a) 



Using the fact that 



4 log det M = Tr flVIM 
dt V 



-1 



one can express the Kolmogorov entropy in the form 



i=0 



lim — - 

T^Qo 2T 







dtlsTii [WUt) + Wlit)] Zi(t) - [Wscit) + Wlit)] Z,{t) 



+ [wut) - wlit)] Z3(t) + [wlit) - w^t)] 1J{t)} 



(42) 
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where we have introduced the operators 



Zi 


= ui'^) 




+ uf ^uf " 




Z2 


= u?) 




+ vf^^vf 




Z3 


= uf 









Taking into account the dynamical equations (pn|) for the blocks of the 
evolution operator, it is easy to see that the operators (pSj) satisfy the differ- 
ential equations 

Zi = £[W,,Zi + ZiW,^, + W,,Z3 + Zl^W,^,-Zi(W,, + W,^,)Zi 

+ zl (W., + Wl) Z3 - Zi (w,, - wl) Z3 - {wl - W,.) Zi] 
Z2 = £[-W,,Z2-Z2Wj,-W,,Zl^-Z3Wf, + Z2(W,, + WjJZ2 

- Z3 (We. + w^J zl^ - Z3 (Wee - wl) Z2 - Z2 {wl - w,.) z^] 

Z3 = e [-W,,Zi - W,eZ3 + Z3W;r, + Z2Wf, - Z3 (We, + Wl) Zi 

+ Z2 (W,e + Wl) Z3 - Z3 (Wee - W^) Z3 - Z2 (W,^, - W,,) Z^] 

(44) 

These equations determine the time evolution of the operators ()43|1 together 
with the initial conditions 

Zi(0) = l, Z2(0) = 0, Z3(0) = 

which can be obtained by substituting in the definitions ()43p the initial con- 
ditions (PT|) for the left blocks of the evolution operator. The system of 
equations pij) can be solved perturbatively by considering solutions of the 
form 

00 

Z,(t) = J]e"Zf)(t) 

n=0 

Substituting these trial solutions into Eq. (jUJ, one obtains the simple result 

Zi(t) = l + o{e) 
Z2{t) = o{e) 

Zsit) = -e [ Wss{r)dT + o{e). 
Jo 

Putting these expressions in Eq. one obtains that the sum of the positive 
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Lyapunov exponents is 



N N N 



J2^^ = 5252^ / {[{Wkn{t)Wkn{t + T)) + {Wnk{t)Wkn{t + T))]cOs{uJn + UJk)T 

+ [{Wkn{t)Wkn{t + r)) - {Wnk{t)Wkn{t + t))] COS K - UJk) t} dr + o{e^) 

(45) 

This formula gives the sum of the positive Lyapunov exponents for the 
dynamical system ()18p. Note that the result has been derived without con- 
sidering any specific form of the coupling matrix Wnk- Therefore it can be 
applied also to cases in which the coupling matrix differs from the form ()19|) . 
In particular, if one considers a symmetric matrix, Wnk = Wkn, the second 
term in the r.h.s. of Eq. (j45j) vanishes and the result reduces to the form 
obtained in ^| for a similar problem. Here, however, the Wnk matrix is not 
fully symmetric, because of the Kronecker delta in Eq. (fT^ which can be 
traced back to the fact that the zero-th Fourier mode has no twin compo- 
nent unlike the other modes (which come in equal pairs [/„ = t/_„). However, 
the second term in Eq. (j45|) can be neglected in the limit of a large number 
of Fourier components, in which case the asymmetry linked to the zero-th 
channel becomes negligible. 

To analyse the physical implications of Eq. ()45p. it is useful to express 
the Kolmogorov entropy in terms of the Fourier components of the random 
potential U. Substituting the matrix elements (fT^ in Eq. one obtains 

Ai = / {Un{t)Un{t + r)) [cos {uJn + UJq) T + COS {uJn " CUq) t] (It 

+ / {Un+k{t)Un+k{t + t)) + {Un-k{t)Un+k{t + t)) COS (u;„ + CUfc) T dr. 

(46) 

Note that we have set the bookkeeping parameter £ = 1, as we will do from 
now on, with the tacit understanding that all results are valid within the 
second-order approximation. Eq. (PUI) represents the Kolmogorov entropy 
for any kind of random potential. We observe that the second sum on the 
r.h.s. of Eq. ()46|) contains a number of terms of order 0{N'^) and therefore 
for iV ^ 1 it is dominant with respect to the first sum which has only 
0{N) terms. Hence if the number of oscillators/modes is large, one can 
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approximate the Kolmogorov entropy (j46|) with 

N N N 



1 r 



(47) 



+ {Un-k{t)Un+k{t + t)) cos {Un + UJk)T dr. 



We remember that the binary correlators which appear in Eq. ()47p are the 
correlators of the Fourier components (jH)) of the random potential. They are 
linked to the binary correlator (jHI) via the double Fourier transform 

{Unmkit + r)) = ^ dy £^ dy' x(r, y - y') cos ^ cos (48) 

If the correlators ()48|) decay sufficiently fast as a function of the difference 
|n — fc| of the indices of the Fourier components of the potential, the second 
correlator in Eq. ()47|) gives a marginal contribution with respect to the first 
one. Therefore one can replace expression (jTfjl with 

V'A, ^ / {Un+k{t)Un+k{t+T)) cos {uJn + UJk)T dr. (49) 

1=0 n=-Nk=-N " 



5.2 Application to specific cases 

We will now apply the general formulae derived in the previous subsection 
to a few particular cases and see how one can recover known specific results 
from the general expressions (j46|) and (j49j) . Let us consider first the case of 
longitudinal disorder, i.e., of a random potential of the form U{x,y) = U{x). 
In this case the only non-vanishing Fourier component (jH)) is the zero-th one 

Unit) = 6noU{t) 

and the binary correlators (PHj) become 

{Un{t)Uk{t + r)) = (T\(r)5„o4o. (50) 

These correlators vanish unless n = k = and therefore there is no doubt 
that they decay fast for increasing values of \n — k\. Hence the approximate 
formula can be applied; inserting expression (jKHjl in Eq. (jinjl one obtains 

N N ^2 roo / 1 \ 

E - E / ^(^) cos 2^;^^ dril - -5no j • (51) 

n=0 n=0 Jo \ / 
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The exact result can be obtained by substituting the correlators (j^UI) in 
Eq. (pUj) : this gives 




cos 2uJnT dr 



(52) 



and by comparing Eq. (jKT|) with Eq. ()52|1 we see that the difference between 
the two expressions is indeed negligible in the large N limit. We remark 
that Eq. ()52|) is in perfect agreement with what one obtains for the sum of 
the positive Lyapunov exponents (jHUjl derived in Sec. |3] From the physical 
point of view, Eq. is a natural consequence of the fact that, as discussed 
in Sec. El in the case of longitudinal disorder the Fourier modes/oscillators 
are decoupled so that the quasi-lD model reduces to a sum of strictly ID 
systems. 

As a special application of formula ()52|) one can consider the case of 
longitudinal white noise. In this case x(r) = S{t) and the Kolmogorov 
entropy (j52j) takes the form 



in agreement with the well-known expression 



for the Lyapunov exponent of a noisy ID oscillator of frequency a;„. 

Eq. ()52|) can be used also to derive the inverse localisation length for a 
ID model, which can be seen as the limit case of a quasi-lD model when 
the transversal dimensions go to zero, L — > 0. In this limit the only elliptic 
Fourier component is the zero-th one and, indeed, setting = in Eq. (j52|) 
one recovers the well-known expression ^ for the inverse localisation length 
in ID models. 

Let us now consider the case of a random potential which is a white noise 
in both the longitudinal and the transversal direction. In this case the binary 
correlator (jH} has the form 
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and the correlators of the Fourier components are 



(f/„(t)f/,(t + r)) 




(53) 



The form of the correlators ()53p implies that the second correlator in the ex- 
pression (jlTj) gives only a marginal contribution, so we can substitute Eq. 
in formula ()49|) . Neglecting terms with 0{N) addends, one obtains 



This equation coincides with the result obtained for a discrete model by 
Hansel and Luciani in jT2| if one identifies the width 2L of the doubled strip 
and the frequencies a;„ in (j54p respectively with the number of channels and 
the square roots of the energies in the model of Hansel and Luciani. 

5.3 A sufficient condition for delocalisation 

The importance of Eq. PU)) lies in the fact that it provides a criterion for the 
onset of a delocalisation transition in quasi- ID models. In fact, the condition 



represents both a necessary and sufficient condition for the suppression of 
orbit instability in the system of parametric oscillators (fTHj) and a sufficient 
condition for the localisation length to diverge. That condition (j55p is not 
necessary for the delocalisation of the electronic states depends on the fact 
that, as discussed in Sec. 12. 6| the localisation length is equal to the smallest 
Lyapunov exponent. Therefore, delocalisation sets in as soon as the minimum 
Lyapunov exponent vanishes, even if larger Lyapunov exponents are non- 
zero. On the contrary, in the case of the dynamical system (jl8p . unless all 
Lyapunov exponents vanish, initially nearby orbits exponentially diverge. 

The importance of expression rests on the fact that it allows one 
to prove that, for specific kinds of long-range correlated disorder, the condi- 
tion is fulfilled over a certain range of the electronic energy and, there- 
fore, a continuum of extended states arises even in quasi-lD models. The 
effect is analogous to the one observed in strictly ID models [21 E]- To 




(54) 



N 




(55) 



i=0 
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see how a delocalisation transition can occur, one can observe that the Kol- 
mogorov entropy (|46p vanishes if the Fourier transforms in the longitudinal 
direction of the binary correlators (jlHjl are zero 



/>oo 

/ (t) (t + r) ) cos ujTdT = (56) 

Jo 

for all values of the indices rii and n2 and for every frequency u in the inter- 
val [0 : 2^/E]. The frequency interval is determined by taking into account 
that the frequencies ()14p vary in the interval [0 : \/E] and that formula pUll 
contains cosines with frequency up to twice the maximum value of the fre- 
quencies ujn- Condition (jSHI), however, is equivalent to the requirement that 
the Fourier transform of the correlator © 



oo 



X{uj^,ujy)= / dx dy x{x,y) cos uj^x cos UyU (57) 
Jo Jo 

should vanish for < tUx < '2.\fE, i.e., 

X(cjx,^s,) = for < cj^ < (58) 

A potential whose binary correlator satisfies this condition can be constructed 
with a slight generalisation of the method discussed in Sec. E] for the longi- 
tudinal disorder case. Starting from a binary correlator (j57j) with arbitrary 
frequency dependence, one can obtain the function 



ii{^,y)= / ^ / -p-'JxiyJx, ^y) exp {-i^x^ - iujyy) . 

J-oo 2tT V 

The required potential can then be constructed via the convolution product 

/■oo rOD 

U{x,y) = a^ dsx I dsy I3{sa:, Sy)ri{s^ + x, Sy + y) 



where rj is a. stochastic process with 

{ri{x,y)) = and {ri{x,y)r]{x',y')) = 6{x - x')6{y - y'). 

Following this procedure, one can obtain a potential which fulfils the 
delocalisation condition ()58|) . For example, one can consider a potential with 
long-range correlations in the longitudinal direction of the form 

1 

x(a;, y) = — (sin z^ix — sin z/2x) 7(?/) (59) 
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where 7(1/) is the correlation function in the transversal direction. The binary 
correlator is the simplest possible generalisation of the correlator 
Its power spectrum has the form 



7(0;^) if Z/i < LlJ^ < Z/2 

otherwise 



By selecting a potential such that ui = 2yE, one is ensured that the delo- 
calisation condition is fulfilled. We are thus led to the conclusion that 
a delocalisation transition can occur if the disorder exhibits long-range cor- 
relations that make the power spectrum vanish in an appropriate frequency 
interval. 

0ns should hasten to add that this conclusion is rigorously true only for 
weak disorder and within the second-order approximation. In ID models 
it has been shown that the delocalisation transition produced by long-range 
correlations of the disorder is a second-order effect and the "extended states" 
are, in fact, electronic states which extend over a spatial range of order 
0(1/(7^) rather than 0(l/cr^) [llj. The same result can be expected for 
quasi- ID models; this does not diminish the practical importance of the 
delocalisation transition analysed here because, for weak disorder, i.e., when 
(T^ — >■ 0, the increase of the spatial range of the electronic wavefunction can be 
huge, being of order 0(1/ct^). For finite samples, therefore, the delocalisation 
can be real and manifest itself in a strong change of the transport properties 
of the disordered sample. 



6 Conclusions 

In this work we have shown how the spatial structure of electronic states in 
a quantum quasi-lD model with weak disorder can be analysed in term of 
the time evolution of a classical system of parametric oscillators with weak 
stochastic couplings. By Fourier-transforming the stationary Schrodinger 
equation for the quasi-lD model in the transversal directions, one obtains a 
set of equations for the Fourier components of the electronic wavefunction 
that can be mapped unto the dynamical equations of a Hamiltonian system 
of coupled oscillators. The spatial behaviour of the Fourier components of the 
wavefunction is thus matched to the time evolution of the oscillators, while 
the disorder in the quasi-lD model manifests itself as noise in its dynamical 
analogue. The specific effect of the noise is to perturb the frequencies of 
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the oscillators and to couple the oscillators among themselves. Both models 
are characterised by a set of characteristic Lyapunov exponents; however, 
whereas in the solid state model the key Lyapunov exponent is the smallest 
one, which is equal to the inverse of the localisation length, in the dynamical 
system the most important exponent is the largest, which defines the mean 
rate of exponential divergence of the orbits. 

The dynamics of the Hamiltonian system can be analysed in full detail 
when the random potential in the quasi- ID model depends only on the longi- 
tudinal coordinate. In this case the oscillators are effectively decoupled and 
it is possible to obtain the whole Lyapunov spectrum of the system. One can 
thus make a complete study of the delocalisation effects produced by specific 
long-range correlation of the disorder. 

The general case of a random potential which depends both on longitu- 
dinal and transversal coordinates is more difficult to handle; nevertheless, in 
it is possible to evaluate perturbatively the sum of the positive Lyapunov 
exponents (or Kolmogorov entropy). Using this result, one can show that 
specific kinds of long-range correlated disorder make all Lyapunov exponents 
vanish within the second-order approximation and therefore produce a delo- 
calisation transition in quasi-lD models like they do in strictly ID models. 
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